On the occurrence and origin of anthropogenic radionuclides found in a fragment of the Chelyabinsk (LL5) meteorite

A piece of the 2013 Chelyabinsk meteorite was investigated for its content of anthropogenic radionuclides. In addition to traces of cesium‐137 that had been previously reported for this particular fragment, we found an unusually high amount of strontium‐90, which indicates that the source of this contamination was the Kyshtym accident (1957). A high Sr‐90/Cs‐137 activity ratio is characteristic for Kyshtym‐derived contaminations. Based on the cesium‐137 content in the soil from the finding site, it is estimated that the fragment was contaminated with soil particles in the milligram range upon impact. Investigation of the soil revealed very unusual ferromagnetic characteristics and an iron‐rich chemical composition. Mössbauer spectroscopy indicated the presence of steel components in this soil, suggesting that the investigated meteorite fragment was found in an industrial dumping site rather than natural soil.     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

Limiting equilibrium method for slope/drilled shaft system

The use of drilled shafts to stabilize an unstable slope has been a widely accepted practice. There are two basic design and analysis issues involved: one is to determine the global factor of safety of the drilled shafts stabilized slope and the other one is to determine the design earth thrust on the drilled shafts for structural design of the shafts. In this paper, a limiting equilibrium method of slices based solution for calculating global factor of safety (FS) of a slope with the presence of a row of drilled shafts is developed. The arching mechanisms due to the presence of the drilled shafts on slope were taken into account by a load transfer factor. The method for calculating the net force applied to the drilled shaft from the soil mass was also developed. The interrelationships among the drilled shaft location on the slope, the load transfer factor, and the global FS of the slope/shaft system were derived utilizing the developed numerical closed‐form solution. An illustrative example is presented to elucidate the use of the solution in optimizing the location of the drilled shafts on slope to achieve the desired global factor of safety of the slope/shaft system. Copyright © 2009 John Wiley & Sons, Ltd.     

Mineral chemical composition and geodynamic significance of peridotites from Nain ophiolite,central Iran

The peridotites from north of the town of Nain in central Iran consist of clinopyroxene-bearing harzburgite and lherzolite with small lenses of dunite and chromitite pods. The lherzolite contains aluminous spinel with a Cr number (Cr^# = Cr/[Cr + Al]) of 0.17. The Cr number of spinels in harzburgite and chromitite is 0.38–0.42 and 0.62, respectively. This shows that the lherzolite and harzburgite resulted from <18% of partial melting of the source materials. The estimated temperature is 1100 ± 200 °C for peridotites, the estimated pressure is <15 ± 2.3 kbar for harzburgites and >16 ± 2.3 kbar for lherzolites and estimated fo₂ is 10^?1±0.5 for peridotites. Discriminant geochemical diagrams based on mineral chemistry of harzburgites indicate a supra-subduction zone (SSZ) to mid-oceanic ridge (MOR) setting for these rocks. On the basis of their Cr^#, the harzburgite and lherzolite spinels are analogous to those from abyssal peridotites and oceanic ophiolites, whereas the chromites in the chromitite (on the basis of Cr^# and boninitic nature of parental melts) resemble those from SSZ ophiolitic sequences. Therefore, the Nain ophiolite complex most likely originated in an oceanic crust related to supra-subduction zone, i.e. back arc basin. Field observations and mineral chemistry of the Nain peridotites, indicating the suture between the central Iran micro-continent (CIM) block and the Sanandaj–Sirjan zone, show that these peridotites mark the site of the Nain–Baft seaway, which opened with a slow rate of ocean-floor spreading behind the Mesozoic arc of the Sanandaj–Sirjan zone as a result of change of Neo Tethyan subduction régime during middle Cretaceous.     

Development in modeling cyclic loading of sands based on kinematic hardening

In this paper, there is presented an elastoplastic constitutive model to predict sandy soils behavior under monotonic and cyclic loadings. This model is based on an existing model (Cambou‐Jafari‐Sidoroff) that takes into account deviatoric and isotropic mechanisms of plasticity. The flow rule used in the deviatoric mechanism is non‐associated and a mixed hardening law controls the evolution of the yield surface. In this research the critical state surface and history surface, which separates the virgin and cyclic states in stress space, are defined. Kinematic hardening modulus and stress–dilatancy law for monotonic and cyclic loadings are effectively modified. With taking hardening modulus as a function of deviatoric and volumetric plastic strain and with defining the history surface and stress reversal, the model has the ability to predict the sandy soils' behavior. All of the model parameters have clear physical meanings and can be determined from usual laboratory tests. In order to validate the model, the results of homogeneous tests on Hostun and Toyoura sands are used. The results of validation show a good capability of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd.     

Historical tsunami in the Makran Subduction Zone off the southern coasts of Iran and Pakistan and results of numerical modeling

Tsunami hazard in the Makran Subduction Zone (MSZ), off the southern coasts of Iran and Pakistan, was studied by numerical modeling of historical tsunami in this region. Although the MSZ triggered the second deadliest tsunami in the Indian Ocean, among those known, the tsunami hazard in this region has yet to be analyzed in detail. This paper reports the results of a risk analysis using five scenario events based on the historic records, and identifies a seismic gap area in western Makran off the southern coast of Iran. This is a possible site for a future large earthquake and tsunami. In addition, we performed numerical modeling to explain some ambiguities in the historical reports. Based on the modeling results, we conclude that either the extreme run-up of 12–15 m assigned for the 1945 Makran tsunami in the historical record was produced by a submarine landslide triggered by the parent earthquake, or that these reports are exaggerated. The other possibility could be the generation of the huge run-up heights by large displacements on splay faults. The results of run-up modeling reveal that a large earthquake and tsunami in the MSZ is capable of producing considerable run-up heights in the far field. Therefore, it is possible that the MSZ was the source of the tsunami encountered by a Portuguese fleet in Dabhul in 1524.     

Spatial distribution and temporal variation of reference evapotranspiration in arid and semi‐arid regions of Iran

Analysis of spatial and temporal variations of reference evapotranspiration (ET_o) is important in arid and semi‐arid regions where water resources are limited. The main aim of this study was to analyse the spatial distribution and the annual, seasonal and monthly trends of the Penman–Monteith ET_o for 21 stations in the arid and semi‐arid regions of Iran. Three statistical tests the Mann‐Kendall, Sen's slope estimator and linear regression were used for the analysis. The analysis revealed that ET_o increased from January to July and deceased from July to December at almost all stations. Additionally, higher annual ET_o values were found in the southeast of the study region and lower values in the northwest of the region. Although the results showed both positive and negative trends in annual ET_o series, ET_o generally increased, significantly so in six (～30%) of the stations. Analysis of the impacts of meteorological variables on the temporal trends of ET_o indicated that the increasing trend of ET_o was most likely due to a significant increase in minimum air temperature, while decreasing trend of ET_o was mainly caused by a significant decrease in wind speed. At the sites where increasing ET_o trends were statistically significant, the rate of increase varied from (+)8·36 mm/year at Mashhad station to (+)31·68 mm/year at Iranshahr station. On average, an increasing trend of (+)4·42 mm/year was obtained for the whole study area during the last four decades. Seasonal and monthly ET_o have also tended to increase at the majority of the stations. The greatest numbers of significant trends were observed in winter on the seasonal time‐scale and in September on the monthly time‐scale. Copyright © 2011 John Wiley & Sons, Ltd.     

List of forthcoming papers

List of forthcoming papers     

Exploring the applications of 3D proximity analysis in a 3D digital cadastre

Increasing population in urban areas and limitations of suitable lands for developing houses and urban infrastructure have led to the vertical development in cities. However, these developments are managed by a cadastral system which is mainly two-dimensional and cannot efficiently represent Rights, Restrictions, and Responsibilities (RRRs) in complex scenarios. In fact, a three-dimensional cadastre is required for efficiently registering and representing RRRs. In this paper, a 3D proximity analysis was proposed and implemented to determine RRRs and associated easement rights in non-topology-based data structures. This method can be used to investigate the surrounding spaces of a subject apartment unit or storage in a high-rise. The performance of the developed method was evaluated in a large complex high-rise in Tehran, Iran. The results confirmed that the proposed method can correctly identify the neighbor spaces in complex scenarios.     

Experimental and 3D Numerical Simulation of Reinforced Shear Joints

The load transfer capacity and failure mechanism of a fully grouted bolt installed across a joint in shear is investigated, both experimentally and numerically, in five types of bolt. The double-shearing testing of bolts were studied in concrete blocks of 20, 40 and 100 MPa strengths, subjected to different pretension loads of 0, 5, 10, 20, 50 and 80 KN, respectively. The parameters examined include: shear resistance, shear displacement, induced strains and stresses during the bolt-bending process, and its ultimate failure across the sheared joint planes. The conclusions drawn from the study were that the strength of the concrete, bolt profile configuration and bolt pretension load played a significant influence on the shear resistance, shear displacement and failure mechanism of the reinforced medium. Experimental and numerical simulations showed that the failure occurs as a result of the induced axial and shear stresses acting between the hinge point distances in the vicinity of the shear joint plane.